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Abstract. The asymptotic behaviour of the superpotential of general SUSY 
transformations for a coupled channel Hamiltonian with different thresholds 
is analyzed. It is shown that asymptotically the superpotential can tend to 
a diagonal matrix with an arbitrary number of positive and negative entries 
depending on the choice of the factorization solution. The transformation 
of the Jost matrix is generalized to "non-conservative" SUSY transformations 
introduced in Sparenberg et al (2006 J. Phys. A: Math. Gen. 39 L639). Applied 
to the zero initial potential the method permits to construct superpartners with 
a nontrivially coupled Jost-matrix. Illustrations are given for two- and three- 
channel cases. 
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1. Introduction 

In the context of quantum scattering inverse problems, not much is known 
about coupled-channel problems with threshold differences, i.e., inelastic problems 
(see paragraph IX. 4 of pQ and references therein). Though well-known single- 
channel methods based on the Marchenko or GePfand-Levitan approaches have been 
generalized to these coupled-channel problems, no sufficiently simple method for 
inverting experimental data has been deduced from these generalizations. This may be 
explained by the fact that in its general form the method requires knowing the whole 
set of scattering data and, in particular, the whole scattering (S) matrix should be 
known, whereas only the open-channel submatrix is accessible from experimental data. 
Another reason is related with the complicated character of the Gel'fand-Levitan- 
Marchenko equation which in this case was solved only numerically [SJ. In the one 
channel case the second difficulty was overcome with the help of the supersymmetric 
quantum mechanics (SUSY QM) approach (see e.g. [31 0]) which is equivalent to 
the Gel'fand-Levitan-Marchenko method when the kernel of the integral equation 
is degenerate |6]. Moreover, when applied to the zero initial potential, these 
transformations lead to exactly-solvable potentials with scattering matrices being 
arbitrary-order rational functions of the wave number and hence providing excellent 
fits of experimental data ,7 . Until recently, however, generalizing such supersymmetric 
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transformations of the zero potential to inelastic coupled-channel problems seemed to 
be impossible as a matter of principle [51 [51 QUI- Thus, experimental data with coupling 
did not seem accessible to this method. 

Fortunately, as it was recently announced [11] , this strong limitation is explained 
by an unnecessary condition imposed on transformation functions used up to now. 
Usual supersymmetric transformations [3l El QI3 can be described in 
terms of transformation operators relating solutions of two Sturm-Liouville problems, 
and hence keeping boundary conditions unchanged; this is why we call them 
"conservative" . Changes in the spectrum due to such transformations may correspond 
only to elements from the kernel of either the transformation operator or its adjoint 
form. This is just the reason why the spectrum of two supersymmetric partners may 
differ only by a finite number of levels. Although the existence of transformations of 
another kind, breaking boundary conditions, is known for some time (see e.g. )12|), 
their use in the single-channel case has been limited to a very specific application 
[13] . where usual transformations happen to be sufficient. In the inelastic coupled- 
channel case, on the contrary, these "non-conservative" transformations seem to be of 
fundamental importance as they lead to exactly-solvable potential models with non- 
trivially coupled S-matrices. We hope that these new transformations may become a 
keystone for an inversion procedure both simple enough for practical realization and 
accurate enough to fit experimental data with a good precision. 

This renewal of interest for supersymmetric transformations of inelastic coupled- 
channel problems implies that a general study of their properties is necessary. In 
particular, the Jost-matrix transformation, from which the transformed S-matrix can 
be deduced, should be known in the general case. In it is shown that the Jost- 
matrix modification depends on the value of the superpotential matrix U{r) at infinity, 
which in turn depends on the asymptotic behaviour of the factorization solution. In 
the single-channel case, U(oo) is a positive (resp. negative) real number when the 
factorization solution increases (resp. decreases) at infinity. In the coupled- channel 
case, U (r) is a matrix and its asymptotic behaviour is much more complicated. The 
main goal of the present work is to study carefully this behaviour. 

Section 2 reviews definitions from multichannel scattering theory. As an 
introduction to Section 3 we first review briefly general properties of multichannel 
supersymmetric transformations and then we analyze the maximal number of arbitrary 
parameters entering into the superpotential and as a result into the transformed 
potential. Next we discuss in a detailed way the superpotential asymptotic behaviour 
proving our main theorem and give a closed expression for the transformed Jost 
solution and Jost function. In Section 4 we derive a general form of the superpotential 
which is a SUSY partner of the zero potential. Section 5 illustrates these findings by 
examples and comparisons with existing results from the literature. Section 6 contains 
conclusions and perspectives. 

2. Multichannel scattering 

Let us first define our notation and briefly recall some notions of scattering theory 
[HI [15] used below. We consider a multichannel radial Schrodinger equation that 
reads in reduced units 

d 2 

Hip(k,r) = k 2 ip(k,r) H = -—r + V (1) 
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where r is the radial coordinate, V is an N x N real symmetric matrix, and ip may 
be either a matrix-valued or a vector-valued solution. By k we denote either a point 
in the space C N , k = {fc 1; . . . , kjy}, k{ E C such that Iraki ^ 0, or a diagonal matrix 
with the non-vanishing entries ki, k — diag(/ci, . . . , fcjv)- The complex wave numbers 
ki are related to the center-of-mass energy E and the channel thresholds Ai, . . . , Ajv, 
which are supposed to be different from each other, by 

*$ = E-A i . (2) 

We do not assume a fixed order of channels since any necessary order can be achieved 
by a suitable permutation of the rows in the Schrodinger equation (Q}. For simplicity 
we will limit ourselves to unequal thresholds here and defer the general study of equal 
and unequal thresholds to a future work. We will also assume potential V to be 
short-ranged at infinity, i.e., there exists an e > such that 

e er \Vij(r)\dr < oo (3) 

o 

where Vij, i,j — 1,...,N are entries of matrix V. Under such assumptions, the 
Schrodinger equation has two Nx N matrix-valued solutions /(±/c, r) (Jost solutions) 
such that 

f(±k,r) — > exp(±i/cr) = diag[exp(±ifcir), . . . , exp(ztifcjvr)]. (4) 

r — >oo 

The columns of these matrices form a basis in the 2iV-dimensional solution space 
of the Schrodinger equation with a given value of E. In general, these solutions 
are complex and satisfy the symmetry property f(k,r) — f*(—k*,r), where asterisk 
denotes complex conjugation; for real energies below all thresholds, k — —k* and the 
Jost solutions are real. 

Next we define the regular solution <p(k, r) and the irregular solution r](k,r) by 
their behaviour at the origin. For the sake of simplicity, we limit ourselves to bounded 
s-wave potentials, in which case these solutions satisfy 

<p(k,0) = (p'{k,0) = I (5) 

rj{k,0) = I r)'(k,0) = (6) 

where prime means derivation with respect to r and / denotes the identity matrix. This 
definition shows that the columns of these matrices also form a basis in the solution 
space of the Schrodinger equation. In terms of the Jost solutions, these solutions read 

<P(k, r) = ^ [f(k, r^Fi-k) - f(-k, r)fc- x F(fc)] (7) 

V(k, r) = i [/(*, r)*- 1 ^-*) - /(-*, lOfc-^fc)] (8) 
where F(k) is the Jost matrix 

F(k) = f T (k,0) (9) 
with T meaning transposition, and matrix G{k) is defined as 

G(fc) = -[/'(fc,0)] T (10) 

Proving ([7]) and ((SJ) requires calculating, both at the origin and at infinity, the 
Wronskian W[<p(k, r), f(k, r)} = (p T (k, r)f'(k, r) — [ip'(k, r)] T f(k, r), which generalizes 
the usual definition of the Wronskian of one-component functions to N channels, 
and the Wronskian W[rj(k,r), f(k,r)]. Equation (JTJ) implies that the value of these 
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Wronskians is independent of r, as well as that of W[f(— k, r), f(k, r)] = 2ik. For 
real energies, both solutions cp and r\ are purely real because they satisfy a system 
of differential equations with real coefficients and real boundary conditions. For 
energies below all thresholds or for energies above all thresholds, this can also be 
directly checked on ((7|) and ((8|), using the symmetry properties F(k) — F*(—k*) and 
G(k) = G*(-k*). 

The Jost matrix defines both scattering and bound states properties. The 
scattering matrix, which is symmetric, reads 

S(k) = k-^Fi-tyF-^k)^ 2 = k 1 ' 2 [F- 1 {k)] T F T {-k)k- 1 ' 2 . (11) 

Bound-state energies, E = E mi m — 1, . . . , M, correspond to zeros of the determinant 
of the Jost function, det-F(iK m ) = 0, such that Re K m > 0, lm.n m ^ — with 
E m — — «m j + A,, i = 1, . . . , N below all thresholds. For potentials satisfying the 
above assumptions, the number M of bound states is finite. 



3. Multichannel SUSY transformations 



According to the multichannel SUSY approach [8l [9], applying the transformation 
operator 

A- = ~ + U(r) (12) 
dr 

to solutions ip(k,r) of ^ leads to solutions ip(k,r) of the new equation 

d 2 

Hip(k,r) = k 2 ip(k,r) H = -— + V(r) (13) 

or 2, 

where V, like V, is supposed to be a real, short-ranged, bounded and symmetric NxN 
matrix and like ip, may be either a matrix-valued or a vector-valued function. 
The matrix-valued function U (usually called superpotential) is expressed in terms 
of a matrix-valued solution of (JTJ) at a fixed value of E = £ below all thresholds 
(this parameter is known as the factorization constant), which we denote a and call 
transformation function or factorization solution. Defining the corresponding wave 
number diagonal matrix n by its positive elements Ki = \/ At — £ , one has 

Ha{r) = -K 2 a{r) (14) 

and 

U{r) =a'(r)o-- 1 (r). (15) 

Then for E ^ £ one has ip = A~ip. The specific form (|12[) . (|T5|) of the transformation 
operator A~ results in the potential V from (| 13[) being of the form 

V ( r ) = V(r) - 2U'(r). (16) 

To have a real and symmetric potential (|16|) we restrict a to be real and such that its 
self-Wronksian vanishes, W(a, a) — 0. 

For E = £ a particular solution of (fT3|) is 

0( K ,r) = (cr^V) H4>(K,r) = -K 2 4>(K,r). (17) 

Other solutions ip corresponding to the same £ may be found as usual from the 
property W(<fi, ip) — —I, which gives 

ip(K,r) = {(T T )- 1 {r) j a T (s)a(s)ds H4>(n,r) = -n 2 ^(n } r). (18) 
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The conventional SUSY transformations have the property that if ^(0) = 
then i/j(0) — {A~i[)) r= Q = which requires some additional limitation on the 
transformation function er(r). Following rejecting this limitation leads to loosing 
this property of the transformation operator A~ . Such transformation operators 
violate the vanishing behaviour of the solution at the origin and we call them "non- 
conservative" . Nevertheless, since A~ transforms solutions of the initial differential 
equation into solutions of the transformed equation the full information about the 
new Hamiltonian H is accessible. In particular its Jost function and 5-matrix can 
be constructed explicitly. Below we will concentrate our attention mainly on non- 
conservative transformations although our main result (Theorem [2] in Section |3 . 2[) is 
valid for the general case. 

3.1. Number of arbitrary parameters in superpotential 

First we discuss the general form of the transformation function a introduced in the 
previous section. We notice that the general vector- valued solution of the Schrodinger 
equation ([I]) contains 2N integration constants and can always be presented as a linear 
combination of 2N fixed linearly independent solutions. In general, transformation 
function a may be composed of N such solutions. Therefore it may contain 2N 2 
arbitrary parameters at most. But as far as the new potential (1161) is concerned there 
is a big redundancy between these parameters. Indeed, because of the specific form of 
the superpotential (| 1 5() a multiplication of the transformation function a on the right 
by a non-singular constant matrix does not affect the superpotential. A minimal set 
of arbitrary parameters is given by 

Theorem 1. Given the initial potential V(r), fixed thresholds and factorization 
energy, the most general transformed potential V(r) is completely determined by the 
value U (0) of the symmetric superpotential matrix at the origin. It is calculated by 
formulas U6\) and i TTSj) where 

<j(r) = 7}(iK, r) + (p(in, r)U(0) (19) 

and contains N(N + l)/2 arbitrary real parameters which are the entries of matrix 
U(0). 

Proof. The general matrix-valued solution of the Schrodinger equation 

a(r) = n(in, r)C a + <p(iK, r)D 1 , (20) 

with real matrices C\ and D\ produces the most general real superpotential (|15[) 
and, hence, potential (fT6)) . Definitions ([5]) and ([6]) then imply that deter(O) = detCi 
vanishes if and only if matrix C\ is singular. In this case cr(0) is not invertible and 
the superpotential and hence the transformed potential V become singular at the 
origin; as stated above, we want to avoid this case here and therefore wc impose the 
condition det C\ ^ 0. It is now clear that in this case matrix C\ does not affect the 
superpotential U as given in (fT5j) since we can multiply ([20)1 by Cf 1 on the right, 
which leaves the superpotential unaffected or, equivalently without loosing generality, 
put Ci = /. The superpotential U and hence the transformed potential V thus only 
depend on the N 2 parameters appearing in D\. 

This simplified writing allows us to express easily the other condition imposed 
on the real transformation function (|20p . namely the symmetry of U and V. As 
mentioned above, this happens when W(a, a) = 0, which gives N(N — l)/2 equations 
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for the elements of matrix D\. The value of this Wronskian being r- independent, (|20p 
can be used to calculate it at the origin with and © which leads to 

W[a(r),a(r)]=D 1 -Df = 0. (21) 

The superpotential and transformed potential are thus symmetric when D\ is chosen 
symmetric. This can also be checked on the value of the superpotential at the origin 
which reads, according to ©, ©, (JTSj) and ((20j) , U(0) = D x . □ 

To calculate the Jost matrix for the transformed potential according to © we need 
to know its Jost solution which is defined by the asymptotic behaviour Usually 
a supersymmetry transformation changes this behaviour. So we have to analyze the 
asymptotic behaviour of the function f(k,r) — A~ f(k,r) which is mainly defined by 
the asymptotics of superpotential U(r). 

3.2. Asymptotic behaviour of superpotential 

According to (|15|). the asymptotic behaviour of the superpotential depends on the 
asymptotic behaviour of the factorization solution. This time, in place of (|2"D)) , we 
choose to write the factorization solution as 

a(r) = /(-uc, r)^ 1 ' 2 ^ + f(in, r) K - 1 / 2 D 2 (22) 

with C 2 and D% being some constant matrices; factor introduced for further 

convenience. We notice that matrices C 2 and D 2 should satisfy the condition 

DjC 2 - ClD 2 = (23) 

following from the symmetry property of the superpotential W[cr(r), c(r)] = 0. 

We show below that the asymptotic behaviour of the superpotential crucially 
depends on the structure of matrix C 2 . In particular, it depends on the rank of 
C 2 and if rankC2 = R < N it depends on an interrelation between the values of 
thresholds and linear dependence between rows of C 2 . This interrelation becomes 
more transparent for a specific order of channels. Therefore before going further we 
will first rearrange the channels taking into account the structure of matrix C 2 . Our 
main aim in this reordering is to collect together both all linearly independent rows of 
matrix C 2 and its linearly independent columns. As it was already noticed changing 
channels corresponds to going to another starting Hamiltonian |T]). But evidently it 
corresponds to the same physical system after the reordering. The permutation of 
columns simultaneously both in C 2 and in D 2 is equivalent to a multiplication on the 
right of the whole factorization solution a by a constant non-singular matrix which 
evidently does not change the superpotential U as given by (fT5|) . 

We rearrange the rows of C% together with the corresponding channels in the 
following way. The first channel with wavenumber k! x and, hence, the first row of 
the reordered matrix (we denote C3) correspond to the largest threshold related to 
a non- vanishing row of C 2 . The second channel with wavenumber k! 2 and, hence, 
the second row of the new matrix C3 correspond to the largest remaining threshold 
related to a row of C 2 linearly independent of the first row of C 2 . At each next step 
i < i?, a new channel with wavenumber k[ and row i of C3 corresponds to the largest 
remaining threshold related to a row of C 2 linearly independent of the previous rows 
of C3. The reordering ends when the first R rows of C3 become linearly independent. 
All remaining rows of C 2 are transferred to C3 without changes. As mentioned above 
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we now permute columns in C3 to have its upper left R x R block non-singular thus 
obtaining matrix C . Matrix D 2 after all these permutations is transformed into D. 
The diagonal wavenumber matrix is written as 

k! 



k" 



(24) 



where k' is the RxR diagonal block after the reordering and k" is the (N—R)x (N — R) 
diagonal block containing the remaining wavcnumbers. Such a structure of C, D and 
wavenumber matrices will be assumed till the end of the paper. 
Now we can formulate our main theorem. 

Theorem 2. When r — > 00, if rank C = R, the superpotential has the asymptotic 
form 

U(r) — * diag(uu) (25) 

uu = +Ki i = l,...,R (26) 

uu = -Ki i = R+l,...,N. (27) 

In order to prove the theorem, we first reduce C and D to canonical forms, 
simplest as far as the superpotential (fT5")) is concerned but reflecting on the one hand 
the singular character of matrix C and on the other hand the non-singularity of the 
whole factorization solution er(r). This is performed in the next lemma. 

Lemma 1. Matrices C and D can be transformed by right multiplication with a non- 
singular square matrix T into the canonical forms 

CT =U 0) DT =( X o -f) < 28 » 

where I denotes the RxR unit matrix in CT and the (N — R) x (N — R) unit matrix 
in DT and Xq is a symmetric matrix, Xq = Xq. Matrix Qq verifies the following 
property. For any i < R and j > R such that the inequality k! 1 < n" holds one has 

q% - (29) 
where q^ are entries of matrix Qq . 

Proof. By construction the RxR upper left block M of C is invertible. Since the 
rank of C is R, the last N — R columns of C are linear combinations of the first R 
ones. This means that there exists an R x [N — R) matrix P such that 

C ( P j ) = (30) 

where / is the (N — R) x (N — R) unit matrix. Similarly, for the rows, there exists an 
(N — R) x R matrix Qo such that 

( Qq -I )C = 0. (31) 

Hence matrix C can be written as 

C = ( Q M Q MP ) (32) 



or 
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The linear dependence of the rows is not modified by the multiplication of C by 
a nonsingular T\ after which row i(< R) contains only one non- vanishing element in 
column i. Property (|29[) now follows from the fact that by construction row j is a 



linear combination of all previous rows i = 1, 2, . . . of C such that k\ > n 
Instead of matrix D 2 the product 



- ( Z Z ) < 35 > 

now appears in (|22[) . According to (|23[) the symmetry property of the superpotential 
is translated into the following conditions: 

£>ii - £>fi = D^Qo - QID21 (36) 
D12 = -QlD 22 . (37) 

Equation (|37p together with form ([34"]) of CTi implies that if D22 in (|35|) is singular, 
function cr(r) given in (|22p becomes singular for all r, a case we would like to avoid 
here so that we necessarily assume -D22 to be invertible. Moreover, using (f3"6")) and 
(|37p we rewrite matrix DTi as 

- - ( x „° ) ( L L ) 

where we denoted Xq = Dn + Finally, we can define T = T\T<z where T\ is 

given by (f3"4"|) and T 2 is the inverted second factor in the right hand side of (j3"5|) which 
keeps unchanged matrix CT\\ the condition Xq = X follows from (|3"6l) . □ 

In the next lemma we establish an important property of matrix-valued function 
Q(r) which will appear in the asymptotic form of the factorization solution. 

Lemma 2. For k', k" and Qo defined as in Lemma{l\ the (N — R) x R matrix 

Q(r) = e K " r Q a e- K ' r (39) 

tends to zero when r tends to infinity. 

Proof. A matrix element of Q is given by 

= &eW-"> (40) 

From Lemma [TJ either the exponential tends to zero or coefficient vanishes. □ 

Now we can prove Theorem [2] 

Proof. According to Qj, the factorization solution has the following asymptotic 
behaviour 

<r(r) - <J as (r)= K - 1 / 2 [c Kr C + e- Kr D]. (41) 

r — >oc 

Recalling that C and D here have canonical forms ([25)1 we get 

i/ 2 f I + X -Q T \ / e KV 



- K 7 I T o e->- » (42) 
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where Q is given by (|39p and matrix 

X(r) = e- K ' r X e- K ' r (43) 
vanishes at infinity. The derivative of (|42|) can be written as 

I-X Q T \ ( e K ' r 



Hence, for U one has the following asymptotic behaviour 



_ 1/2/ Q T \ ( i + x -Q T \ 1 1/2 



From Lemma [2] and (|43|) one obtains for r — ► oo 

MO - K f f °A (46) 



x " J , 

which concludes the proof. □ 

Corollary 1. If the thresholds are ordered such that K\ > «2 > ■ ■ • > and 
matrix C has rank R the number of arbitrary parameters in the superpotential is 
R(R+l)/2 + R(N-R). 

Proof. For the given order of thresholds all may be chosen different of zero 
according to Lemma [TJ Moreover, the parameters enter in the superpotential only 
through matrices X and Qo- □ 

3.3. lost-matrix transformation 

After the asymptotics of the superpotential is found we are able to calculate both the 
Jost solution and the Jost matrix for the transformed potential. 

Theorem 3. The Jost matrix F(k) of the transformed potential V reads, in terms of 
the Jost matrix F(k) and the function G{k) H10\) of the initial potential, 

F{k) = [U(oo) - ifcp 1 [F(k)U(0) + G(k)} . (47) 
Proof. According to Theorem [2] the Jost solution of the transformed potential reads 

f(k,r)=A-f(k,r)[U( X >)-ik}- 1 (48) 
as seen with (Tl"2"|) and (UJ). The theorem then follows from definitions © and (|10[) . □ 
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4. SUSY partners of V{r) = 

The zero initial potential is important since in this case compact analytic expressions 
are possible both for the transformed potential and for its Jost function. The initial 
Jost solution in this case is simply the exponential f(k,r) = cxp(ifcr) and the initial 
Jost function is the identity matrix, F(k) = I. Hence, G(k) = — ifc. The regular 
solution has the form: ip(k,r) — sin(fcr)fc -1 and the irregular solution is written as 
rj{k,r) = cos(fcr). The factorization solution respecting Theorem [1] has the form 

a(r) = cosh(Kr) + sinh(Kr)K _1 t/(0) (49) 
= \ c Kr [I + K-'-UiO)] + \ e— [I - K~ l U(0)] . (50) 
Another parametrization corresponds to Lemma [1] 

a(r) = n- 1 ' 2 [e Kr C + e - Kr D] (51) 
with an appropriate choice of matrices C and D. 

Proposition 1. Let I be the R x R identity matrix, matrices Xq and Qo be chosen 
according to Lemma{J\ Q(r) be defined by iS9\) and X(r) by |^3| ). // the parameters 
are such that detY(r) =/= W G [0,oo) where 

Y{r)=I + X{r) + Q T {r)Q{r) (52) 

then the potential V(r) = —2U'(r) with 

is a (non-conservative) SUSY partner ofV(r) = and has the Jost solution 

f(k, r) = [U(r) - ik] e lkr \U{oo) - ifc] -1 (54) 

and the Jost function 

F{k) = [U{oo)-ikY 1 [U{{))-\k] (55) 

where U(oo) may be found from H25\\27\ ). 

Proof. First we notice that the function a(r) = <r as (r) given in (|4ip is just function 
(|5ip and, hence, it can be taken as transformation function to produce a SUSY partner 
of V(r) = 0. For C and D in block forms (|28| the superpotential is given in (|45"1) . 
From here after some algebra one gets (|53|) . Expressions (|54|) and (|55|) for the Jost 
solution and the Jost function correspond to (|4"5|) and (|4"T]) for the zero initial potential, 



respectively. □ 
Corollary 2. For Xq = the superpotential L53\) can be written as 

U = k- 2k 1 /2qT (Q r Ql) - 1 Q r «Va ( 56 ) 

= - K + 2K 1 ' 2 Q c (Q T c Q c y 1 Q^k 1 ' 2 (57) 



where Q r and Q c are row and column block matrices written in terms of Q t39\) as 
Qr = (Q, —I) where I is the (N — R) X (N — R) identity matrix, = (J, Q T ) where 
I is the R x R identity matrix. 
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Proof. Using the property 

Q[I + Q T QY 1 =[I + QQ T Y X Q (58) 
one obtains from (l53l) 



QQ T / \ I + QQ T ) (59) 

= K " 2K (, -Q / ){ 7 + 0Q T J (60) 

Equations (fBT)]) and (f5"T)) are nothing but compact forms of ([50)1 and (fSU)) 
respectively. □ 

For two particular cases corresponding to rank C — N — 1 and rank C = 1 either 
(|56[) or (|57p takes a particularly simple form since in one case Q r is a row and in the 
other case Q c is a column. The explicit expressions are given in the following corollary: 

Corollary 3. Let channels be ordered such that K\ > K2 > ■ ■ ■ > Let also 

Q be a row, Q — (q%, . . . , g^v-i), q,; = q® exp(ftjv — Ki)r where qf = for any i such 
that kn > Ki and arbitrary otherwise. Then the superpotential has the following block 
form: 

TT u- 2 ,„l/2 ( Q T Q ~Q T V.l/2 fa-,\ 

U - K -TT0Q? K { -Q i ) K (61) 

where Q T Q is an (N — 1) x (N — 1) matrix with entries qtqj, i,j = l,...,N—l. 

Let us now order only the first channel such that Ki > max(K2, ■ ■ • , kn) and Q T 
be a row Q T = (q2, . . . ,qN), q% — <Zi exp(/Cj — K\)r. Then the superpotential may be 
written as 

^-' t Hk"'1J«')• , " (62) 

where QQ T is an (N — 1) x (N — 1) matrix with entries qiqj, i,j = 2,...,N. 

Proof. The statement follows from Lemma [T] and Corollary [2] For the first part of the 
statement rank C — N — 1 whereas for the second part rank (7=1. □ 

Another simplification occurs for rank C = N and a particular choice of matrix 

X . 

Corollary 4. Let Xq be a column of N arbitrary real parameters and X = 
exp(— nr)Xo. Then the superpotential reads 

U = K- YT ^ K ^XX T n^. (63) 

Proof. Choosing rank C = N (meaning that Qo = 0) and Xq = XqXq , we have 
Y = I + XX T . Using the property (XX T f = (X T X)XX T one gets Y- 1 = 
I - (1 + X T X)- 1 XX T . The statement follows now from (J53J). □ 

5. Examples 



Let us now illustrate the theorems and Proposition[T]we have just established by some 
exactly-solvable examples, supersymmetric partners of V(r) = 0. 
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5.1. The 2x2 model with rank C = 2 (Cox potential) 

Let us start from a two-channel problem. 

According to Theorem Q] we choose the transformation function as given in ([50)) 
with the maximal number of arbitrary parameters included in 

U(0) = ( a l f V (64) 



[3 a 2 

Taking (|5U|) into account, the condition rank C = 2 reads 

(Ki+ ai ){K 2 + a 2 )-0 2 ^0, (65) 

Because of the simple character of the transformation function we easily find the 
superpotential according to (|15p . Its off diagonal elements have the form: ii 12 = 
u 2i = (3/ deter. For the first diagonal element one obtains 

cosher) 

tin = — [ax costu^r) + Kx smh(K 1 r) 

det a 

[(aia2 — /3 2 ) cosh(Kir) + Ki«2 sinh(Kir)] . (66) 



Here 



K2 det a 

012 . O-X 

deter = — cosh(Kir) sinh(K2r) H cosher) sinh(/cir) 

K 2 Kx 
(XxOL2 — (3 2 

+ cosh(Kir) cosh^^) H sinh(Kir) sinh(K2r). (67) 

K\K2 

The element m 22 is obtained from (|66p by the replacement kx ^2 and ai <-» «2- 
According to (JT6J) the transformed potential is simply twice the derivative of these 
expressions with the opposite sign. 

Since rank C = 2, Theorem [2] states that no negative entries in the asymptotic 
form of the superpotential can appear so that U(oo) = k = diag (ki, /ta) which can 
also be checked by a direct calculation. Applying now Theorem [3] we find the Jost 
matrix for this potential 

F(k) = ^JlL (68) 

\ K 2 — ifc 2 K2— 1^2 / 

We notice that up to a change of parameters we obtain one of the two-channel 
Bargmann potentials previously found by other means [16) . 

The same superpotential may be rewritten with another parametrization. 
According to Proposition [T] for rankC = N the superpotential reads U = —n + 
2k 1 / 2 (I + X)^ 1 n 1 / 2 with X given by ([43]) the elements of which are iy — exp(— n^r — 
Kjr)xfj, i,j — 1,2. Its more explicit form is 

1+X22 -XX2 
(1 + Xxx)(l + X22) ~ xl 2 " V _X 12 1 + XX1 

After the change of the parameters 

*Sx = det[U ' 0) + K] \P - (<* - + "»)] (70) 

^ 2 = ZTTTTT^ , n [/3 2 - + «i)(a 2 - « 2 )] (71) 



U--K+ - k 1/2 ^ 22 12 U 1 / 2 ("691 



det[t/(0) + k] 

-2^ v /k7k 2 
det[£/(0) + k] 



*? 2 = ~frf^, (72) 
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one recovers the previous result. 

Two explicit examples of this 2x2 model are given in figures Q] and [5] Part (a) of 
each figure shows the potential while part (b) shows the eigenphase shifts 6%, 62 and 
the mixing parameter e [TS]. For these examples, we have chosen the thresholds 

Ai = A = 10, A 2 = 0, (73) 

which implies that the factorization energy £ and wave number {fti,Ka} are related 
by 

n\ = &-£> k% = -£. (74) 
In both cases, the value of the superpotential at the origin, (|64p. is chosen as 

^°)=(o1 -2)' ( 75 ) 

while the factorization wave number is K2 — 3 and K2 = 2.2 in figures [T] and 
respectively. Comparison of figures [1] (a) and fj] (a) shows that the choice of 
factorization energy strongly modifies the potential. In contrast, figures [T] (b) and 
[2] (b) show that this choice does not strongly affect the scattering matrix, which 
always displays a Feshbach resonance at an energy of about 6.3 and a large negative 
slope of 82 at zero energy and of 5\ above threshold. The eigenphase shifts and mixing 
parameter are only a bit smaller in figure [2] (b) than in figure [1] (b). 
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5.2. The 2x2 model with rank C = 1 

In the previous N = 2 example, matrix C has rank 2. Let us now consider the 
case, not allowed in [16j . where rank C is 1. Our supersymmetric formalism, on 
the contrary, is also valid in this case. Both matrices Qo and Xq defined in Lemma 
[1] are numbers, Qo = QOi Xq = xq, so that Q(r) = q(r) = qo exp(K 2 r — Kir) and 
X(r) = x(r) — xq exp(— 2kit). Choosing K\ > k 2 and xq > —1 — <7q we get from (|53[) 



with 



l-x~q 2 l + x-q 2 2< ? ^/k 1 k 2 
ai = ;Ki a 2 = — t K2 o=- -. (77) 

It is not difficult to check the condition rank C = 1 (cf. (p5]) ) 

(Ki+ai)( K2 + a 2 )-/3 2 = (78) 

where ai j2 = 01,2(0) and (3 = fe(0). We notice that our approach permits us to 
calculate the superpotential by the same formulas (155)) and ([57]) as in the previous 
section but now the parameters are not independent anymore, they should satisfy 
condition (|78[). With the identification xq = and 4kik 2 //5 2 = (go + l/<7o) 2 this leads 
to an exactly-solvable model used in [11] to construct an analytical model for the 
Feshbach-resonance phenomenon. 

Theorem [T] implies that U(oo) — diag(Ki,— k 2 ) which is clearly seen from (|76p . 
The value U (0) is given by the same formula (|76[) with the replacement x — » xq and 
q —t qQ. Therefore applying Proposition [T] one gets the Jost matrix 

( ai-ifci \ 

which differs from l|68p because of the different asymptotic form of the superpotential 
at infinity. 

Let us now construct an explicit example that illustrates the strong impact of 
U{oo) on the scattering matrix. In figure [3l a potential is constructed with the same 
parameters as in figures [1] and [21 except for the factorization wave number which is 
now chosen in order to satisfy condition ([75)1 . k 2 w 2.194 675 . . . When k 2 is smaller 
than this limit, the potential becomes singular, a case we want to avoid here. Figures[T] 
(a),H](a) and[3](a) show that when /t 2 approaches this limit value from above, the well 
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in potential V22 goes to infinity and finally disappears when k 2 actually reaches the 
limit value. As long as n 2 > 2.194 675 . . ., this well movement in V22 has practically 
no impact on the scattering matrix: the potentials are nearly phase equivalent with 
each other. In contrast, when the well disappears, a strong change of behaviour is 
observed, as seen in figure [3] (b): while the Feshbach resonance still keeps the same 
energy and width, the slope of 82 at zero energy becomes small and positive. Above 
threshold, 5 2 has now a value very close to 5\ of figures [1] (b) and [2] (b) , with a large 
negative slope. Though 82 is continuous at threshold, it now has a strong cusp effect. 
The mixing parameter, which was close to zero, now gets close to —tt/2, while 8\ gets 
very small. The rank of C, and hence !7(oo), thus have a strong qualitative effect 
on the scattering matrix, which displays very different behaviours in both cases; for 
practical applications, e.g. for fitting actual scattering data with such potentials, both 
behaviours might be of interest. 

Let us finally remark that the behaviour U{po) = diag(Ki,— K2) agrees with 
the result obtained in [5] for N — 2 in the context of a bound-state removal by 
supersymmetric quantum mechanics. There, the factorization matrix solution is made 
of one increasing and one decreasing vector solution at infinity, which corresponds to 
rank C = 1. In [8], this result is even generalized to N channels and a diagonal matrix 
for U(oo) is found, with all positive elements but one, corresponding to the channel 
with the lowest threshold. However, the general situation is more complicated, as 
described by our Theorem [2] there may be more than one negative element in the 
asymptotic form of the superpotential and the lowest channel may correspond to a 
positive asymptotic entry of the superpotential. This possibility is illustrated by our 
last example. 



5.3. The 3x3 model with rank C = 2 

Here we choose K\ > K3 > K2, rank C = 2 so that matrix Qq is the row Qq — (qo, 0). 
Such a choice of Qq reflects the fact that the third row of C is proportional to its 
first row (see (f29|)). The adopted order of thresholds corresponds to k' — diag (k\, K2) 
and k" — K3. For simplicity only (equal) off-diagonal entries of 2 x 2 matrix Xo from 
Proposition [T] are chosen different of zero which we denote xq. Then applying (|53|) we 
obtain 

2 ( 1 -x q \ 

U = -k+- k 1/2 -x l + q 2 -xq n 1/2 (80) 

deter I 9 / 

\ q -xq q- ) 

where deter = 1 + q 2 — x 2 , x — xoexp(— Kir — K2f), q = goexp^r — Kir). The 
condition det a 7^ Vr € [0, 00) is satisfied if q 2 . > Xq — 1. It is clearly seen here that at 
r — > 00 both q and x vanish and the asymptotic form of U, C/ as = diag(Ki, K2, —Ka), 
agrees with Theorem [21 



6. Conclusion and perspectives 

In conclusion, supersymmetric quantum mechanics is a promising tool to pragmatically 
solve the inelastic coupled-channel inverse problem, thanks to its ability to construct 
sophisticated exactly-solvable models. The present work is only a starting point to 
the general study of supersymmetric transformations in the multichannel case: by 
focusing on the asymptotic behaviour of the superpotential, it reveals the richness of 
possible behaviours, as compared with the single-channel case. 
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We have made a rigorous study of this asymptotic behaviour for a general 
supersymmetric transformation. Our theorem generalizes a result found in the 
literature [5J [TO] : there may be several negative elements in the asymptotic form of 
the superpotential and the lowest threshold may produce either a positive (as shown 
explicitly by our example 1 5. 3 1) or a negative entry. 

As a byproduct of our proof, we have been able to construct exactly-solvable 
potentials that are supersymmetric partners of the zero potential. Essential 
simplifications occur when each element of the factorization matrix solution is either 
an increasing or a decreasing exponential (not a linear combination of both types). 
Explicit examples have been given for two and three channels. 

Future work should focus on both the behaviour of the superpotential at the origin 
and at infinity in the general case: arbitrary number of channels, with equal and/or 
different thresholds, and arbitrary rank for the matrices multiplying the regular and 
singular solutions in the factorization matrix solution. On the other hand, iterations 
of coupled-channel supersymmetric transformations should be studied, which should 
eventually lead to a satisfactory solution of the coupled-channel inverse problem, with 
and without threshold difference. 
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